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ABSTRACT 


In this thesis, a new VSAT (Very Small Aperture Terminal) communication 
system is developed using Direct Sequence Code Division Multiple Access (DS-CDMA) 
for multiple maritime mobile users in the footprint of national communication satellites. 
The Forward Error Correction (EEC) is implemented by applying convolution encoding 
with soft decision decoding. The worst-case scenario is always considered by placing 
VSAT system and different types of jammers on the footprint where the minimum signal- 
to-noise ratio is possible. Using this assumption, the performance of the system is 
analyzed for different convolution code rates, for a different number of users and for the 
different jammer powers. The Walsh codes are used to establish an orthogonal cover 
between CDMA channels in a VSAT. Additionally, spread spectrum is included by PN 
sequences to ensure as much orthogonal coverage as possible between the VSATs. Also, 
that spreading is the key to minimizing these hostile jammer effects. After analyzing 
VSAT system for possible scenario elements, optimum system parameters are introduced 
for military and civilian applications. 
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I. INTRODUCTION 


A. BACKGROUND 

Mankind fought with each other long before recorded history. As the centuries 
passed, new long-range weapons appeared on battlefields and on seas. Long-range 
weapons gave opponents ability to hide and to make surprise attacks, but at the same time 
they created the need to locate the enemy. This need eventually led us to the invention of 
radar. While World War II was raging in Europe and later in other places as well, the 
techniques and capabilities of radar usage was upgraded. 

The first use of electromagnetic action in history began in World War I and in the 
Russian-Japanese War, long before the “Electronic Warfare” concept was shaped. On 
May 31, 1916, Admiral of the British Fleet, Sir Henry Jackson, used coastal radio 
direction finders to observe movements of the German fleet and, on the basis of this 
information, positioned the British fleet so that it could oppose the German forces 
successfully. Earlier use of electromagnetic action is claimed where the Russian torpedo 
boat Gromky and cmiser Igumrad jammed Japanese radio communications during the 
Russian-Japanese War at the Battle of Tsushima in May 1905 [1]. 

As a principle of radar, the reflected energy that is returned to the radar from a 
target indicates the presence of a target, but also transmitting such electromagnetic 
energy reveals presence of a radar to the target. These types of intentional and 
electromagnetic actions and counter actions occured many times in many campaigns of 
WWH. Meanwhile, all types of electromagnetic actions began to be called “Electronic 
Warfare.” 

The “electronic warfare” concept was shaped during WWH. During the Cold 
War, the potential effect of Russian EW systems over U.S. weapon systems strengthened 
the importance of EW. Admiral Thomas H. Moorer, Chairman of the U.S. Joint Chiefs of 
Staff in 1970-1974, stated, “If there is a WW IE, the winner will be the side that can best 
control and manage the electromagnetic spectrum” [1]. 
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Today we must use the electromagnetic spectrum well and protect our weapon 
systems from enemy electromagnetic actions to achieve success. A reliable 
communication is also essential for both sides. Since this communication system uses 

the electromagnetic spectrum, it will also be targeted by the enemy’s electronic warfare 
capabilities. Until we find another way to communicate other than radio waves, 
communication signals will intentionally be jammed by the enemy or be interfered 

unintentionally by friends. At this point, there is an urgent need for a reliable 
communication system, which can be used in all weather conditions, is endurable against 
enemy jammer activity, and is difficult for an enemy to intercept. 

B. OBJECTIVE 

In this thesis, we will design a satellite communication system with a “Very Small 
Aperture Terminal” (VSAT) that can easily be carried by a naval vessel. This system will 
serve the naval platforms in the footprint of national communication satellites to carry all 

of the voice traffic as well as the data. Data traffic will essentially consist of position 

information of friends and verified enemies, which will be detected by sensors aboard 
ships or helicopters. Continuous transmission of this valuable information from all 
platforms in the tactical arena would be important to winning the battle. In some cases, 
interruption in transmission cannot be tolerated even if E mission Control (EMCON) 

measures are on progress. Due to the directivity advantage of satellite terminal antennas, 
a minimum electromagnetic signature would be revealed to the enemy outside of the 
narrow main beam, so continuous transmission could be established in many cases. To be 
able to establish a continuous picture of the tactical arena, some technical requirements 
must be achieved, such as bit rate and bandwidth for desired signal-to-noise plus 
interference ratio that can be supported by the platform. 

This system includes many tools such as Direct-Sequence Code Division 
Multiplexing Access (DS-CDMA) and convolutional encoding with soft-decision 
decoding. Different ratio of contributions from spectrum spreading and convolutional 
encoding for a fixed spreading factor enables us to have a flexible jammer resistant 
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system against different capabilities of the enemy jammer. Significantly, we will analyze 
the effect of enemy jammer types on our system and define the capability of this system. 


C. RELATED WORK 

In this thesis, we mainly studied the performance analysis of spread spectrum 
satellite communication systems under different jamming conditions. A detailed similar 
related work was done in [3], studying the ‘Modeling and Analyzing of Cellular CDMA 
Forward Channel.” Also in reference [15], Ugur Yigit has studied the performance 
analysis of a spread spectrum VS AT system for different coding techniques. Chapter in 
of this thesis and of reference [15] were written jointly by the author and U. Yigit. 


D. THESIS OUTLINE 

The main goal in Chapter II is calculating the minimum signal-to-noise ratio that 
our system can support. Initially, we define the term VSAT and other orbital mechanics 
parameters. Then we decide on the proposed VSAT system requirements. By using the 
physical sizes, distances and power limitations, a link budget analysis will be introduced. 
The worst-case scenario is considered in this analysis assuming a friendly platform at the 
edge of the satellite footprint and the enemy jammer at the best place to minimize our 
system’s performance. The link budget analysis presents the minimum and maximum 
signal-to-noise ratios that can be established at the satellite receiver. The difference 
between the minimum and maximum signal-to-noise ratios will be calculated by the link 
margin analysis. 

In Chapter HI, we examine the algebraic properties of the Walsh functions and PN 
sequences that will help to determine how to minimize jammer effects and friendly 
interference on our system. We present the extended orthogonality that ensures the 
orthogonality between every channel in a VSAT and introduce the PN sequences used to 
spread out our frequency spectrum, reducing the power spectral density and minimizing 
the jammer effects. Afterwards, we introduce the convolutional forward-error-correction 
coding, which allows the system perform better with a low signal-to-noise ratio. 
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The performance analyses will be introduced in Chapter IV. The performance 
analyses are considered both for tone jammer and pulse jammer. First of all, we develop a 
signal-to-noise plus interference ratio and the probability of bit error for our DS-CDMA 
satellite system by using Gaussian approximation. For both jammer cases, co-channel 
interference, inter-VSAT interference, Additive White Gaussian Noise and Jammer 
Interference were considered as variances of our random variable Y . In our uplink 
model, we introduce Forward Error Correction (FEC) in the form of convolutional 
encoding with soft-decision and develop an upper bound on the probability of bit error 

[3]. 

Then, we use a simulation software to analyze the probability of bit error for the 
minimum signal-to- noise ratio, which will be defined in Chapter n. At the end of 
Chapter IV, optimal performance parameters will be introduced. 

We will summarize the conclusions in Chapter V. 
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II. VERY SMALL APERTURE TERMINAL 


In this chapter, basic satellite communication and VSAT concepts will be 
introduced. Some necessary parameters in satellite communication and orbital mechanics 
will be assessed to find the maximum signal-to-noise ratio for the modeled VSAT 
system. Then rain attenuation will be calculated to decide on the link margin for the 
minimum uplink signal-to-noise ratio. Minimum and maximum SNR values will be used 
in Chapter IV to decide on whether reliable communication can be achieved or not. The 
decision given on reliability was based on the probability of bit error for a given signal- 
to-noise ratio range. We will accept that the probability of bit error should be below 10 6 
within minimum and maximum SNR values to establish reliable communication. 

A. DEFINITION OF VSAT 

Since Marconi made the first radio transmission across the Atlantic in 1901, 
communication has been one of the fastest developing branches of electrical engineering. 
If military requirements had not been so vital when radio was invented, we might never 
have had a chance to use modem communication devices today. Unquestioningly, 
military needs have led the developments of communication techniques and devices at 
the times of WWI, WWII. and Cold War. 

While communication technology meets the demands of the military and civilian 
community, each advancement in these technologies encourages the military and civilian 
sectors to demand even more sophisticated technology. Consequently, new technology 
continually begets more and more progressive applications. Thus, today, modem military 
commands, on the field, in the sea, and in the air require that vast packets of data be 
transferred rapidly (through high bit rates), and further require large bandwidths, 
mobility, power efficiency, secrecy and durability under friendly or hostile interference. 

A dream of communicating with the farthest point on earth or even across the 
solar system came within reach after the first communication satellite was set in orbit. 
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VS AT systems are merely a link in this chain. A VS AT system is only an example, which 
can answer all those requirements listed above. 

New developments in military command structure, which now use information 
and intelligence more than ever before, exchange enormous amounts of data between 
different commands and units. In peacetime, establishing several communication media 
that enable such operations is possible. However in a crisis situation, establishing an 
effective communication system in the field, either in one’s own territory or in a foreign 
country is troublesome to maintain. An affordable solution that can resolve most of these 
problems is the VS AT. 

The term VSAT stands for Very Small Aperture Terminal. VSATs are small 
ground terminals, which can communicate with a satellite either as a receive-only 
terminal or a transmit-receive terminal. VSATs provide cost-effective solutions for the 
growing telecommunication needs throughout the world [10]. Today’s satellites are more 
powerful, enabling the use of smaller and less-expensive antennas on the ground. Also, 
the developments in silicon technology include most of the necessary VSAT functions, 
which makes VSATs more effective. 

Although there are many definitions for a VSAT, the European 
Telecommunication Standard Institute has a comprehensive one, which mainly describes 
a VSAT as part of a satellite transmit-receive system that has an aperture size smaller 
than 3.8 m 2 [11]. 

The main VSAT deployment configurations are the star and mesh configurations. 
Each system has advantages and disadvantages. For military applications, the mesh 
structure has more advantages than the star configuration. The central hub in a star 
configuration would be an important and primary target for an attack, which would cause 
an overall catastrophe for communication if the central hub were destroyed or were made 
to malfunction. There is no vital VSAT in the mesh configuration. Therefore, even if the 
high traffic hubs all malfunctioned, the rest of the system would be able to carry the 
entire communication load, not through the hubs, but directly from VSAT to VSAT via 
the communication satellite. 
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B. LINK GEOMETRY 


Today's VSATs are commonly but not necessarily used with geostationary 
communication satellites. A geostationary satellite always appears to be at the same 
altitude and azimuth for a particular location on earth. Therefore such a satellite rotates 
around the earth with a period of one sideral day and with zero inclination with respect to 
the earth’s equator and stays above the same geographical location as the earth rotates at 
the same velocity as well. The gravitational force upon a body decreases as the distance 
from the earth increases. The centrifugal force increases upon a body as the velocity of an 
orbiting body increases. At a certain distance - a point at which the satellite is precisely 
located - the gravitational pull toward the earth’s center and the centrifugal force away 
from the earth become equal for a circular orbit. This distance is known as geostationary 
radius, a gso . The geostationary radius is 42,164 km regardless of the size or of the mass 

of the satellite or orbiting body. The derivation of this distance was given in reference 
[8]. The geostationary height is found by subtracting the earth’s radius from a gso . 

Ko =a S so ~ a e 

= 42164-6378 (2.1) 

= 35786 km 

In the link budget analysis, we will use the slant range from the VS AT to the 
satellite to calculate the path loss. The geostationary height and the slant range (distance 
from the satellite to user) are only equal to each other at the sub-satellite point. The sub- 
satellite point is defined as "the point on earth vertically under the satellite" in reference 
[8]. The sub-satellite point is the projection point of the satellite on the equator for 
geostationary satellites. Because of this fact, the geostationary satellite positions are 
defined as the longitude of the sub-satellite point (() vv ). 

A worst-case scenario was considered in both the link budget analysis and in the 
performance analysis. We design the model for the farthest user in the footprint. As 
shown in the EIRP map (Figure 1), the footprint of Turksat 2A reaches far beyond the 
Mediterranean to the Atlantic Ocean. We will set the borders of the coverage area as the 
48 dBW transmitter EIRP (where GIT = 4 dB / K at the receiver) contour for the optimal 
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usage. We created the optimal contour for the following criteria: the optimum contour 
must encircle a considerable amount of operational waters for the Tur ki sh Navy with a 
maximum ( GIT) k ratio. The 48 dBW contour encircles the Black Sea, the Agean Sea 

and 90 % of the Mediterranean Sea. In our model the user will be on the 48 dBW contour 
and be positioned as far west as possible. This position was determined to be 

38.00° N -0 °E for the position of the earth station (X E =38 and <|) E = 0). 



Figure 1. Turksat 2-A West Fixed Beam Footprint and EIRP Map (From [7]) 


The slant range derivation for 38.00° N -0 °E position is introduced as below [8]: 


d sr = yja E + a 2 gso -2a E a gso cos (b) 

(2.2) 

b = ai'ccos (cos((|) £ -§ ss )cos(^ £ )) 

(2.3) 
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Substituting for b in (2.2), 


d sr ^<3 £ + a gso 2a E a gso cos((f) £ )cos(A, £ ) 

= ^6378 +42164 2 -2x6378x42164x cos(0-42)cos(38) 
= 38247.15 km. 


As seen in the result, the slant range for our particular user is almost 2,500 km 
longer than the geostationary orbit altitude. This 2,500 km of difference affects the path 
loss considerably. 
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C. VSAT SYSTEM REQUIREMENTS 

1. Link Budget Analysis 

In a possible conflict, such as in the Mediterranean, different types of platforms 
will be used in the scenario. Each platform will have different communication needs, 
such as different bandwidths or different bit rates. For instance, a task force would need 
the sum of all the task group bit-rate requirements in order to relay necessary reports 
between the task groups and the command base. So this system must be able to serve 
different platforms, which have different properties. 

Older ground terminals used Frequency Modulation (FM) for communication 
with a satellite. Today, with the advances in digital signal processing, much more 
effective schemes can be used. These new methods enable increased data rate, increased 
reliability and efficiency. One particular method is Code Division Multiple Access 
(CDMA). CDMA allows multiple users to share the same bandwidth and still respond to 
the varying bandwidth requirements. It also allows integrating the users with small 
bandwidth requirements into a high data-rate satellite network [12]. 

Our primary objective is to configure a VSAT (mesh structure) based satellite 
communication system, which can maintain ship-to-ship, ship-to-shore and shore-to-ship 
communication via satellite. Because this system is proposed to be on ship, there will be 
no need to make a power efficient system. But still a tradeoff regarding the power to be 
used in this system is inevitable. The more power used, the more the bit error rate 
decreases. On the other hand, the more power used, the more possible interception by an 
enemy becomes. 

Finally, the system should achieve a data rate of at least 6.312 Mbps and must be 
able to maintain a connection with the satellite under adverse weather conditions and/or 
jamming. Using spread spectrum techniques enhanced with proper error correction 
coding algorithms will make our VSAT communication system resistant to friendly 
interference and hostile jamming. The required signal-to-noise ratio calculation for the 
uplink without jammer contribution is necessary for the further performance analysis. 
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A link analysis is given for the uplink according to specifications for Turksat 2 A. 
We will not need the downlink calculation for the system because the performance 
analysis will be based on the scenario that the communication satellite is the target to be 
jammed. 


Satellite Name 

Turksat 2-A 

Origin 

Turkey 

Mission 

Telecommunications 

Location 

42 E (Geostationary) 

EIRP (Fixed West Beam) 

53 dBW (Max) 

48 dBW (Optimum coverage) 

Satellite G/T 

9 dB/K (Max) 

4 dB/K (Optimum coverage) 

Uplink Frequency 

17300-18100 MHz 

Downlink Frequency 

11700-12500 MHz 


Table 1. Specifications of Turksat 2-A (From [7]) 


Proposed VSAT specifications are given in Table 2. 


Required Bit Rate 

6.312 Mbps (T2) 

VSAT Power 

150 W 

VSAT G/T max 

20 dB/K 


Table 2. Proposed VSAT Specifications (From [17]) 


The minimum VSAT Antenna Gain for uplink was found by taking the minimum 
carrier frequency into account: 


G = 


4t iAJ 2 471 X7tx0.8 2 x 0.55x(l7300xl0 6 )- 
c 2 (3xl0 8 ) 2 


= 46211.71 = 46.64r/5 (2.4) 
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Effective Isotropic Radiated Power: 


EIRP = P I G, 

= 150x46211.71 
= 6931757.42 
= 68.4 MB 


Channel (Path Loss): 


L,(path loss) 


4tt fD 


V c J 

^471 xl7300xl0 6 x 38247 xio 3 a2 
3xl0 8 


( 2 . 6 ) 


= 7.68x 10 20 = 208.85<r/5 


Hie maximum SNR for the uplink at the maximum slant range is; 

= EIRP(dB) + — (dB) - L c (dB) u - L (dB) -R h (dB ) +228.6 (2.7) 

where L s : System Losses 

System losses include Receiver Feeder Loss (RFL), Antenna Misalignment Loss 
(AML), Atmospheric Absorption (AA) and Polarization Loss (PL). Assuming the worst 
scenario, we will accept that RFL = 2 dB, AML = 1 dB, AA = 2 dB and PL = 1 dB. 

Therefore the total system losses will be added to 

L s = RFL + AML+ AA + PL 

= 2 + 1 + 2 + 1 ( 2 . 8 ) 
= 6 dB 

The maximum Signal-to-Noise Ratio therefore is 

= 68.41+4- 208.85 -6 -101og6.312xl0 6 + 228.6 =18.16dB (2.9) 
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2 . 


Rain Loss 


Hie most important reason to look for a suitable rain loss is to ensure that the 
communication system will be in service without being interrupted due to inevitable 
events that do not always exist. For example, rain is a factor that must be considered in 
link margin analysis. The effects of rain are accepted as the major problem in satellite 
communication at frequencies above 10 GHz [9]. Unlike receiver feeder loss, antenna 
misalignment loss or polarization loss, rain attenuation is not included in the system 
losses. All the components in the system’s losses can be minimized, even completely 
e lim inated. Conversely rain attenuation is unexpected and inevitable; therefore, it can 
neither be e lim inated nor minimized. Applying a link margin is basically designing the 
system better than it really needs to be under normal conditions. The designer must be 
very precise at this stage. Considering uplink, any unnecessary use of sources can cost 
much and reduce benefits. Considering down link, an unnecessarily established link 
margin can reduce the satellite’s life span and cause the satellite to drop out off service 
earlier. 



Figure 2. Crane Global Model Rain Climate Regions for Europe (From [9]) 
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Figure 2 illustrates the different rain precipitations for different locations in 
Europe and in the Mediterranean. As we set up our model over maritime conditions, only 
the Mediterranean and Black Sea will be considered. 

Most of the Mediterranean is covered by D i region and the Black Sea is covered 
by F region. Since attenuation is severe in the Mediterranean, I) l rate distribution values 
will be used in rain loss analysis. 


Rain Rate 

Exceeded 

Percent of 
Year 

Point Rain Rate Distribution Values (mni/hr) per Rain 
Climate Regions 

D, 

d 2 

d 3 

F 

0.001 

90 

102 

127 

66 

0.002 

72 

86 

107 

51 

0.005 

50 

64 

81 

34 

0.01 

37 

49 

63 

23 

0.02 

27 

35 

48 

14 

0.05 

16 

22 

31 

8.0 

0.1 

11 

15 

22 

5.5 

0.2 

7.5 

9.5 

14 

3.8 

0.5 

4.0 

5.2 

7.0 

2.4 

1.0 

2.2 

3.0 

4.0 

1.7 

2.0 

1.3 

1.8 

2.5 

1.1 


Table 3. Point Rain Rate Distributions for the Rain Climate Regions of the 

Crane Global Model. From ([9]) 


Rain attenuation is defined in terms of exceeded percentage of the time within a 
year and the rain rate corresponding to that exceeded period [9]. For instance, rain 
attenuation is exceeded for 0.1 percent of the time in a year, for a rain rate of 11 mm/hr 
precipitation. In this case, the rain rate is defined as R 0 , where p = 0.1. 


a = a R h p clB I km 

(2.10) 

\ = a K L s r „ dB 

(2.11) 
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where a and b depend on frequency and polarization. The values of “a” and “b” are 
published in reference [9]. Considering uplink frequency (17.3 GHz) values of “a” and 
“b” are a h =0.0367, a v =0.0335, b h = 1.154 and b v =1.128 for 15 GHz. Subscript “h” 
and “v” are for horizontal and vertical respectively. L, is the physical distance in the 
transmission path that radio waves are attenuated. This distance is a function of the 
antenna look angle and the thickness of the raining layer. 


L _ K~K 


sin (El) 


( 2 . 12 ) 


where 

h r =rainheight 

h 0 = VSATheight(sea level) 

El = AntennaLookAngleorSatelliteElevation 

The rain height can be found in reference [8] as a function of the latitude and the 
rain rate. In Figure 4.4 of reference [8], the rain height was found 3.5 kn for 0.1 % rain 
rate and latitudes between 30 to 40 degrees. 


7 . 


El = arccos 


= arccos 


-sin (arccos (cos((|) £ -(f> s;s )cos(X £ ))) 


d 

42164 . 


38247 


sin (arccos (cos(-42)cos38))) 


(2.13) 


26.61° 


L _ El—JL 

sin (El) 

3.5 

sin(26.67) 

= 1.8 km 

Finally, the last term in (2.11) is the reduction factor r . The reduction factor is 
calculated with the equation extracted from reference [9]. 
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90 


90 + 4L i cos(£7) 

__90_ 

~~ 90+4x7.8cos(26.67) 
=0.76 


(2.14) 


Hie common “a” and “b” values can be found as, 

a h +a v 0.0367+0.0335 

a , = —-- =-= 0.0351 

2 2 

b _ a h b h + a A _ 0-0423 + 0.0378 _ 1 m 
c 2 a c 0.0702 

Therefore the total rain attenuation is 

A p =0.0351x1 1 114 I x7.8x0.76 km 
= 3.2dB 


(2.15) 

(2.16) 


Of course the rain attenuation is not the only factor in the link margin analysis. 
Water vapor and other gas and particle compounds in the atmosphere effect the link 
margin. Ippolito has published a detailed study of radiowave propagation in satellite 
communication in reference [9]. 

The maximum SNR was introduced in (2.9). The minimum SNR is found by 
subtracting the rain attenuation from the maximum SNR. 



= EIRP(dB) + —(clB)-L c (dB) u - L s (dB)-R h (dB) +228.6- L m 

min s 


= 18.16-3.2 dB 
~ 15 dB 


16 



In any unexpected, extraordinary weather conditions, performance becomes worse 
due to the exceeded link margin tolerances. In such a case, the output power of VSAT 
can be increased to compensate for the unexpected conditions. For instance, doubling 
output power increases SNR 3dB. In the presence of a pulse jammer ( E b / J 0 = 10 clB ) for 
the code rate of 1/4, a 3dB increase in the SNR would decrease the probability of bit 
error from 1 x 10 4 to 7 x 10 1 . 

D. SUMMARY 

In this chapter, we define the term VSAT in antenna size terms. According to the 
given definition, the satellite communication system must have an antenna smaller than 

3.8 m 2 to be called a VSAT. In other words, the antenna diameter should be less than 2.2 
m. In the model that we present, the antenna size was 1.6 m in diameter. Naval ships can 
easily carry and maintain such an antenna. 

Additionally we introduced a scenario inside the footprint of the Turksat 2-A 
communication satellite. A friendly ship was positioned at the edge of the coverage area. 
Pessimistically, the lowest GIT for the satellite receiver and the farthest slant distance 
were assumed. System losses were exaggerated and rain attenuation was considered. 
Thereafter, the signal-to-noise ratio and the rain loss were calculated for given 
parameters. These calculation results will be used in the performance analysis in Chapter 

IV. For a reliable voice communication, the probability of bit error of 1CT 3 can be 
sufficient in civilian use. But in military communication systems, this probability is 

desired to be lower than 10 6 . 

In the next chapter the proposed DS-CDMA model will be introduced by 
comparing the similarities and differences between cellular and satellite communication. 
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HI. VS AT UPLINK MODEL 


In this chapter we will introduce the Walsh functions and PN sequences that will 
be used in proposed VS AT system together. We will achieve CDMA and rejection of co¬ 
channel interference between channels in the same VSAT by using Walsh functions. 
Likewise, PN sequences will minimize the inter-VSAT interference between the same 
channel sets that are used in different VSATs. The implementation of this system is 
presented in Figure 3. 


b x {t) -► 


“► to the transmitter 


C\ (t ) Wj(t) cos(2t i f c t) 


Figure 3. Spreading of Signal with PN Sequences c(t) and Walsh Function w(t) 

A. WALSH FUNCTIONS COMBINED WITH PN SEQUENCES 

Walsh functions mean everything for digital DS-CDMA communication systems, 
such as satellite or cellular communication. In cellular communication, each unique 
Walsh code in a Walsh sequence provides orthogonal cover on the forward traffic 
channel within each cell to e lim inate intracell interference. Each mobile user in the cell 
has a unique Walsh function assigned, which encodes the traffic coming from the base 
station to the mobile handset. By applying his unique Walsh function to the traffic 
coming from the base station, the mobile user despreads and, in effect, decodes only the 
traffic that is intended for a specified user. When his unique code is applied to the 
remaining intracell traffic in the channel, the orthogonality between user’s code and other 
user codes zeros out their interference. A simplified representation of this process is 
depicted in Figure 4 [3]. 
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Figure 4. Direct Sequence Spread Spectrum (From [3]) 

On the other hand, Walsh codes are not the only way to spread the signal. 
Pseudorandom noise sequences are used to spread the signal as well. Spread 
Pseudorandom noise is a non-orthogonal function so the coverage supplied by the PN 
sequence would not be as successful as coverage with a Walsh code. Tlie intracell traffic 
carrying information for other users would remain spread across the channel's frequency 
band causing interference for the other users in adjacent cells due to imperfect 
orthogonality between the cell’s PN sequences. How this interference affects the 
performance of the communication system will be expressed in detail in the next chapter. 
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f 

Figure 5. Effect of PN Sequence on Transmit Spectrum 


f 


We have introduced an existing example of Walsh codes and PN sequences 
application, which is cellular communication. Considering the VS AT system, Walsh 
coding and PN sequences can be applied to satellite communications with some 
differences in modeling. In cellular communications, every user communicates with his 
or her own base station in the same cell with different Walsh codes and the same PN 
sequence. Similarly, in satellite communications, every VS AT communicates with a 
single communication satellite with one PN sequence but channels separate from each 
other by Walsh codes. In other words, in a VSAT model every terminal represents cells, 
which are separated by different PN sequences. Every channel on the single VSAT will 
be representing users, which are separated by different Walsh Codes. 

1. Walsh Codes 

Walsh functions have many desirable properties, which we will examine in this 
section. We will look at their autocorrelation functions, and orthogonality, which are 
important to the communications engineer. 

a. Properties of Walsh Function 

Walsh functions are generated from generator matrices, Hadamard 

matrices, Rademacher functions and Walsh binary index. Orthogonal Walsh functions are 
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defined in order of N as W N ={wj(t);te (0,T), j = 0,1,..., iV — l}, consisting of N = 2 k 
elements that are functions of time and that have the following properties [13]. 

• Wj(t) takes on the values {+1,-1} except at a finite number of 
points of discontinuity; where it is defined to be zero. 

• Wj(t) = 1 for y' = 0 

• Wj(t) has precisely j sign changes in the interval (0.T). 

T r [0, if j + k 

• f wj($)w k ($)dt = \ . (3.1) 

J o L T, if j = k 

(Orthogonality Property) 

In Figure 6, a set of Walsh functions of order eight has been illustrated for 

W a = { w 0 00 ,w,( t), w 2 (t ), w 3 (f), w 4 (0, w 5 (0, w 6 (i ), w 7 (0} 

In this case, the interval (0, T) has been broken into N = 8 pieces, each 
section being T c =T / N=T/ 8 time units long. 
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Figure 6. Set of Walsh Function Order Eight (From [3]) 


The most significant property of the Walsh functions as applied to the DS- 
CDMA communication systems is that of orthogonality between two different Walsh 
sequences. We will use this property frequently in our performance analysis. Details 
about analysis will be introduced in the next section. The orthogonality property implies 
that the integrating the product of any two different Walsh sequences over a period is 
always zero. An explanation of the fourth property is depicted below. 
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w 3 (0 = [l 1 0 0 1 1 0 0] 
w 6 (f) = [ 1 0 10 0 10 1] 
w 3 (t).w 6 (t) = [0 110 10 0 1] 

Values of "-1" in Figure 6, indicate logic zero in this example where logic 
one is indicated with "+1." The product of w 3 (t).w 6 (t) is calculated by the "exclusive or " 

(XOR) function. The XOR function outputs logic one for the different inputs and logic 
zero for same inputs. Since there are equal numbers of ones and zeros in the product 
matrix of w 3 (t).w 6 (t), integration over a period is always zero. 

b. Extended Orthogonality 

Orthogonality may not always be maintained with three Walsh sequences 
from a set. The result of the integration of two different Walsh sequences may not be 
orthogonal to the third Walsh sequence. 

J 0 w i W w , 0) w k (0 =j 0 w k OK (0 = T (3.2) 

An example of the situation given in Equation (3.2) is given below for the 
Walsh function set in Figure 6. Notice that the integration of w 6 {t).w 1 (f).w l {t) over a 
period is "T" 


w 6 (t ) =[ 10100101 ] 
w 7 (t) =[10101010] 
w 6 (t).w 7 (t) =[0 0 0 0 1 1 1 1] 

Wj(f) =[l 1 1 1 0 0 0 0] 

w 6 (t).w 7 (t).w 1 (t)=[ 1 1111111] 

One way around this problem is to dim inate some of the Walsh sequences 
from the complete set of function W N and force Equation (3.3) to be valid for every 

Walsh function in the reduced set W N , which we use to define extended orthogonality. 

J 0 w > w j ^ w k < 3 ) =J 0 W I ( ? ) K0 = 0 (3.3) 
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Besides the advantage of establishing orthogonality for every channel on VSAT, reducing 
the number of Walsh sequences in a set also reduces the number of channels [3]. In 
extended orthogonality, we will be limiting the number of channels per VSAT. In the real 
world, some VSAT terminals may require more channels than others do in overall 
network. In such cases, a greater order of Walsh functions must be used to achieve 
sufficient number of extended orthogonality Walsh sequences. Then, the system can 
support a sufficient number of channels. The primary central hubs in the star 
configuration or the auxiliary hubs definitely need more complex Walsh functions than a 
remote terminal. Methods for generating a maximum set of extended orthogonality Walsh 
functions are given in reference [3]. 

c. Autocorrelation of Walsh Function 

In this section, we will examine the autocorrelation functions for the 
Walsh functions and average autocorrelation for a set of Walsh functions. In order to 
explore the autocorrelation function of Walsh functions, we will extend our definition 
and consider each Walsh function to be periodic with a period of T , which is consistent 
with their use in practice. Each channel on VSAT has a specific Walsh sequence. A 
satellite onboard processor encodes the data by applying the channel's entire Walsh 

sequence. The Walsh functions are periodic signals with a period of T . Accordingly, the 
Walsh function and the corresponding autocorrelation function is also periodic with 
period T . 

We define the normalized autocorrelation function, a,(n), for any 

periodic Walsh function w l (t) <= W N by 

1 T 

a Xu) =— [w.(t)w.(t-u)dt. (3.4) 

T J 
i o 

Although Equation (3.4) is used for continuous signals, we can consider 
Walsh sequences as if they were discrete signal. This allows us to generate the 

autocorrelation of the Walsh functions easily. In spread spectrum systems, the chip 

duration, T , is defined to be the bit duration divided by the number of chips per bit. 
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(3.5) 



We will consider T as the time increments in discrete Walsh sequences to build the 

autocorrelation functions. The Walsh function presented in Figure 6 is now presented as 
the autocorrelation function in Figure 7. 


1 

OV,(") o 



Figure 7. Set of Autocorrelation Function for W 8 (From [3]) 
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The superposition of all the autocorrelation sequences in Figure 7 divided 


by N gives the "Average Autocorrelation Function" as: 

4v(*0 = 77 £“/(“) (3-6) 

where a,(w)was defined in (3.4). The average autocorrelation function of 
is generated from (3.6). 



Figure 8 . Average Autocorrelation Function, \(u) for the set W N (From [3]) 


We find that the form of A, a (u) is similar to the form of the autocorrelation 
function for a random binary signal. A random binary signal is generated from an 
infinitely long binary random sequence in which the bits are independent and identically 
distributed random variables. The resulting norm aliz ed autocorrelation function for the 
random binary signal is defined after [13] by 


P(«) = * 


\u\ 

\u\ 

T~’ 



0 otherwise 


(3.7) 


2. PN Sequences 

A pseudorandom noise (PN) sequence is a sequence of binary numbers, which 
appears to be random, but in fact is perfectly deterministic. The PN spreading on the 
uplink provides near-orthogonality and, hence, min im al interference between signals 
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from each channels sets from different VSATs. This allows reuse of the band of 


frequencies available, a major advantage of the CDMA. The Direct Sequence Spread 
Spectrum (DSSS) uses a secondary modulation, faster than the information bit rate, to 
spread the frequency domain content over a larger band. The spreading process not only 
assures the system to be more jammer resistant, but also ensures it not to be detected by 
enemy interceptors. The Low Probability Intercept (LPI) systems use spread spectrum as 
well. In DSSS, each data bit is modulated by a Pseudonoise (PN) sequence that 
accomplishes the spectral spreading. The PN sequence consists of random- lik e plus and 
minus ones, which are called “chips.” Each data bit is modulated with at least 11 to 16 
chips. Therefore modulated data bits seem like PN sequences in time domain after being 
modulated. 

a. Properties of PN Sequences 

The ratio of data bit intervals to the chip durations are known as 
processing gain. The higher the processing gain, the better the autocorrelation properties, 
and hence, the better the ability to reject narrowband interference. In Figure 9-a PN 
sequence modulates the data bits, 10 10, (Figure 9-b), resulting in the modulated data 
bits in Figure 9-c. 


■ i i i i * i t i i ” 
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Figure 9. Modulation with PN Sequence 

In this example, the processing gain is calculated as 13 because there are 
thirteen chips per bit. The demodulation process is established in the satellite receiver by 
applying the exact PN sequence at the correct time (Figure 10). Hie timing process will 
be introduced in the following section. 
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Figure 10. Correct Demodulation with PN sequence 


Data bits (10 10) were recovered after the applying the PN sequence to 
the modulated data bits. Figure 9 has illustrated the time domain characteristics of a DS 
Spread Spectrum signal. The frequency domain representation is given Figure 11. The 
over modulation of the data signal leads to a lower power spectral density covering a 
larger frequency band. Suppose the spread signal is transmitted in the presence of a 
narrowband jammer, the despreading operation at the receiver will take the wide band 
spread spectrum signal and collapse it back to the original data bandwidth. The receiver 
will also act on the narrowband jammer so that its spectrum is spread and causes much 
lower interference to the despread signal. This is known as jamming resistance or the 
natural interference immunity of spread spectrum signals. 

We can see that DS Spread Spectrum is bandwidth inefficient in that it 
uses N chips to transmit a single bit of information. Without spreading the spectrum, we 
could have transmitted N bits in the same bandwidth. This inefficiency is the tradeoff to 
achieve interference rejection, or the ability to have reliable communications even in the 
presence of an interfering signal, such as a jammer. It also reduces the power spectral 
density of the transmitted signal so that its transmission causes less interference in other 
systems operating at the same time on the same frequency band. 
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Recovered Data Signal 






Figure 11. Direct Sequence Spread Spectrum 
b. Autocorrelation of PN Sequence 


PN sequences have good autocorrelation properties to allow the receiver to 
recover bit timing and to synchronize the receiver with the uplink. In satellite 
communication, we expect the users transmit spontaneously without any synchronization 
between platforms. In such a case, a satellite receiver must apply the Walsh Sequence and 
PN sequence in the exact time not to allow phase difference between the incoming signal 
and the demodulation signals. The satellite receiver decides the perfect timing by using 
the autocorrelation of the incoming signal. The autocorrelation function outputs the peak 
value at the moment of the correct timing. Figure 9 has illustrated the recovering data 
signal by applying the correct PN sequence without any phase difference. If the receiver 
did not apply the correct timing to demodulate the incoming signal, the data would not be 
recovered. 

«•) 

b) 

c) 

d) 

Figure 12. Incorrect Demodulation of Data Signal 
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An illustration of such a case is given in Figure 12. In part a) is the PN 
sequence that is used in Figure 9. Instead of multiplying the PN sequence (Figure 12-a) 
and modulated data bits (Figure 12-c), we multiply the shifted version of the PN 
sequence (Figure 12-b). The shifted version of the PN sequence was derived by shifting 
the PN sequence one chip left (early). The result was plotted in Figure 12-d. The decision 
mechanism after the PN demodulation is different from the demodulation with the Walsh 
functions. Due to perfect orthogonality, the integration of the demodulated signal over a 
period has given either zero or one as the decision statistic. However, the PN sequences 
do not support full orthogonal coverage. Therefore, we must use a decision mechanism 
other than the one in the Walsh function. We check the after demodulation-integration 
results with a predefined threshold, which is zero, to decide on the data bit. Although the 
transmitted and modulated data sequence was (10 10), the decision b made as (0 1 0 1) 
after the integration in the example given in Figure 9. 

B. CONVOLUTIONAL CODING AND ENCODING [15] 

The digital communications itself is prone to errors. The information transferred 
between the satellite and earth stations is affected by external factors and depending on 
the signal strength, errors are introduced. The amount of error that can be tolerated 
depends on the application. For voice communications, bit error rates (BER) up to 
10 'can be acceptable, but data communications require BER of at least 10 6 . Using 
higher signal power can increase the performance, but this itself alone is not enough to 
correct all the errors. Also, the high power may not be available everywhere and the non- 
linearities in amplifiers li mit the output power. Because of these reasons, error correction 
mechanisms are used. 

The principle of error correction coding is adding redundant bits to the 
information bits and using these redundant bits to detect and correct the errors at the 
receiver. 
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Redundant Bits 



Figure 13. Convolutional Encoding 
The code rate is defined as: 

_ k 
n 


(3.8) 


n : Number of encoded bits 


k : Information Bits 

Therefore, in a coded system, there are (n-k) redundant bits. The coding is 

usually introduced before the modulation in the transmitter, and decoding is done after 
the demodulator at the receiver. In 1948, Shannon demonstrated that by proper encoding, 
errors in the received sequence could be reduced to a desired level without sacrificing the 
rate of information transfer. Provided that for an AWGN channel, Shannon’s capacity 
formula is 


f p ^ 


( S \ 

1 + 

N„Bj 

= 5 log. 

1 + — 


C = B- log 2 

C : Channel Capacity (bits per second) 

B : Transmission Bandwidth (Hz) 

P : Received Signal Power (W) 

N 0 : Single-Sided Noise Power Density (W/Hz) 


C . 

— = log, 
B 


\ + R ^ 


V 


BN n 


(3.9) 
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where P = E b R b and C /B is the bandwidth efficiency. 

There are two methods for coding: Convolutional Coding and Block Coding. The 
block codes are linear and calculate the output frame by depending on the current input 
frame only, so each block is coded independendy of the others and it has no memory. The 
convolutional codes store the memory of previous input frames and use this to encode the 
current input frame. Convolutional codes have finite memory and they are linear also. 



Information 

Figure 14. Encoding and Decoding Process 

The errors can be corrected once they are detected. The detection of errors is 
achieved by the redundant bits. Once the error is detected, it can be corrected in several 
ways. The most common methods are “Error Concealment,” “Automatic Repeat 
Request” (ARQ) and “Forward Error Correction” (EEC). 

EEC tries to recover the original information from the received data. A simple 
example of a (5, 2) block code can be used to demonstrate the FEC: 


Information 

Code 

00 

00000 

01 

01011 

10 

10101 

11 

11110 

Table 4. Encoding 
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A look-up table or a logic construction can be used to implement the values. In 
the Table 4, at least three bits would have to change in any code sequence to produce 
another code sequence. Therefore, this code has a minimum distance of three. 

When a code sequence is received, it is verified for its correctness. If there are 
errors in this received sequence, the code compares it with other sequences and tries to 
correct it. It is clear that the code can easily correct 1-bit errors. For example, if the 
transmitted sequence was 01011 and it is received as 01001, the received sequence will 
have a difference of one bit from 01011, two bits from 00000, and three bits from 10101 
and 11110. Thus, it can be corrected as 01011. On the other hand, 2 bits of error will be 
closer to another sequence, and it will be miscorrected by the decoder. In some cases, the 
detector can recognize that more than one error has occurred and declare that there are 
uncorrectable errors. The FEC can be used with other correction mechanisms in 
conjunction, resulting in even lower error rates. 

The cost of error correction coding is the decrease in the information bandwidth. 
If we represent the uncoded bit rate as R h and the coded bit rate as R coded , the code rate 
is: 


r is always less than 1. For constant carrier power, the bit energy is inversely 
proportional to bit rate. Therefore, 

E 

coded 

for BPSK modulation, the BER is 


eyncoded 


Q 


' N n 


For a coded bit stream 


eroded 
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v N » , 
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This shows that P efoded is larger than the P emcoded , i.e. the probability of bit error 

with coding is worse than without coding. However, this is the error rate at the input of 
the decoder. When the demodulated bit streams are fed into a decoder, some of these 
errors will be corrected and there will be a coding gain. 

A convolutional code is generated by passing the information bits through a finite 
state shift register. The shift register consists of N ( k -bit) stages and n linear algebraic 
function generators, as shown in Figure 15. The input data is shifted into and along the 
register, k bits at a time. The number of output bits for each k bit input data sequence is 
n bits. The number N , called the “constraint length,” indicates the number of input data 
bits the current output is dependent upon. The constraint length determines the 
complexity and the power of the code. Convolutional codes can be described by their 
generator matrices, tree diagrams, trellis diagrams or state diagrams. 


◄- N stages-► 



Figure 15. Encoded Sequence (N Bits) 


Decoding of convolutional codes is a more difficult problem than encoding. The 
function of a convolutional decoder is estimating the encoded input information using a 
method that results in the minimum possible number of errors. U nlik e a block code, a 
convolutional code is a finite state machine. Therefore, the output decoder is a 
“maximum lik elihood estimator” and optimum decoding is done by searching through the 
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trellis for the most probable sequence. Depending on whether hard decision or soft 
decision decoding used, either Hamming or Euclidian metric is used, respectively. 


Rappaport has given the following definition in [16] for the Viterbi algorithm, 
which is valid for both hard and soft decision decoding: 

Let the trellis node corresponding to state S f at time i be denoted 5 •,. Each node 

in the trellis is to be assigned a value V(S j ( ) based on a metric. The node values are 
computed in the following manner: 

1. Set V(S 0 0 ) = 0 and i =1 

2. At time i, compute the partial path metrics for all paths entering each node. 

3. Set V(Sj ( ) equal to the smallest partial path metric entering the node 
corresponding to state S f at time i. Ties can be broken by previous node choosing a path 
randomly. The non-surviving branches are deleted from the trellis. In this manner, a 
group of minimum paths is created from S Q 0 . 

4. If i <L + m , where L is the number of input code segments (k bits for each 

segment) and m is the length of the longest shift register in the encoder, let i = i +1 and 
go back to step 2. 

Once all the node values have been computed, start at state S 0 , time 
i = L + m, and follow the surviving branches backward through the trellis. 

The resultant path is the decoded output for the input stream. 


C. SUMMARY 

In this chapter, we introduced three types of coverage over our satellite uplink 
transmission. The Walsh functions supply full orthogonal coverage, which is capable of 
eliminating all inter-VS AT interference in the uplink. The only problem with the Walsh 
coverage is the possibility of having a non-zero result for the more than two sequence’s 
product integration. The extended orthogonality method compensates for this weak side 
in the Walsh coverage with an expense of a less number of usable Walsh sequences. 

Modulating the uplink signal with a PN sequence makes our system spread in 
spectrum and ensures the system can resist interference. As an effect of spectral 
spreading, low power spectral density prevents the VSAT system from being intercepted 
by enemy ESM capabilities. Even though the PN sequence does not support full 
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orthogonal coverage, it prevents the interference between different VSATs, which are 
using the same set of Walsh functions. However, completely eliminating interference 
between the channels using the same Walsh sequence in different VSATs is not possible. 

The coverage with the PN sequences also establishes a processing gain, also 
called the PN sequence spreading factor. Multiplying the PN sequence spreading factor 
with the processing gain from the convolutional coding results in the overall spreading 
factor of our proposed system. We have previously stated that the overall spreading factor 
of the uplink was 256. 

Digital communication systems are subject to noise due to natural reasons. Error 
correction mechanisms can reduce the required signal-to-noise ratio for a desired BER. 
There has always been a tradeoff between the coding gain and the bit-rate. System 
designers must chose the most appropriate code rate and bit-rate pair according to their 
requirements. We introduce a variable code rate, which can be compensated with 
processing gain from PN spreading. 

We present the connections between Spread Spectrum Systems, Walsh Sequences 
and Orthogonality, PN sequences and Error Correction Coding in this chapter. When all 
of them are combined, a robust and high performance system can be designed. The next 
chapter analyses a system that uses all these components. 
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IV. PERFORMANCE ANALYSIS OF DS-CDMA SATELLITE 

UPLINK 


This chapter examines two types of jammer effects over a VSAT system. These 
are pulse jammer and tone jammer. We will introduce the probability of bit error solution 
for a spread spectrum VSAT system using convolutional coding. According to these 
solutions, performances of a VSAT system against different jammer types will be 
compared. 

A. PERFORMANCE ANALYSIS WITH TONE JAMMER 

A tone jammer transmits an unmodulated carrier with power J somewhere in the 
spread spectrum signal bandwidth. The one-sided power spectral density of tone jammer 
is shown in Figure 16. The single-tone jammer is important because the jamming signal 
is easy to generate and rather effective against direct sequence spread spectmm systems. 
Analyzing this jammer with coherent SS systems shows that the jammer tone should be 
placed at the center of the SS signal bandwidth to achieve maximum effectiveness [2]. 

As we indicate the similarities and differences of cellular and satellite 
communication, the total number of M VSATs carrying K channels are placed in 
footprint of a communication satellite (Figure 17). The number of channels per VSAT is 
not a fixed number for each VSAT system. Some of the terminals can have less or more 
channels than the majority according to tactical needs. Besides the VSATs in the 
footprint, a hostile jammer will be in the middle of the footprint to achieve a maximum 
effect on our communication system. In this worst-case scenario, we will assume that our 
friendly VSATs are placed at the edge of the footprint where the sensitivity of the 
receiver is least. 
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Figure 16. Typical Jammer One-sided Power Spectral Densities (From [2]) 




Figure 17. Multi-channel Allocation 
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We will attempt to recover the information bits /?, (t ) from a specified channel 

(k =1) on a VSAT (i =1). The received signal from all channels and VSATs will be 
despread, demodulated and finally integrated over the bit to form a decision statistic, Y . 
A single tone at the carrier frequency is used for coherent detection and for extracting the 
demodulated signal, y 2 (t) [3]. 

r(t) = S 0 (t ) + n(t) +q (t) + J (t) (4.1) 

where, 

S 0 (t) : Information + Co-channel Interference, 

S 0 =I 0 + y 0 (4-2) 

n (t) : Additive White Gauissian Noise 

q (t) : Inter-VSAT Interference 

J (, t ) : Jammer Interference 

A (t) y 2 (t) 

'■« — *-0 —* 0 —* 0 —► 

AAA 


c(t ) w t (0 COS(27t f c t) 



Figure 18. Despreading and Coherent Detection of Incoming Signal from VSATs 

1. The Despread Signal, y x (t) 

The despread signal is obtained by applying the appropriate PN sequence, which 
is synchronized with the transmitting VSAT. Synchronization is made by using the 
autocorrelation function of the PN function, as described in Chapter HI. The next step in 
the despreading is applying the Walsh sequence in order to remove the orthogonal cover 
as well as to remove spreading. 
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y 1 (t) = r(t)c( t) m{( f) 


y 1 (t)=S 0 (t)c(t)w l (t)+n(t)(( t) w x (t) +c,(t)c(t) w x (t) + J (t) c( t) w;( t) (4.3) 

Hie equation above can be simplified as follows: The subscript “1” indicates the 
despreading. 


Yi(0 - W +yi(0+n l (t)+q l (t) + J l (t) 


K -1 

I l (t)+y l (t) = E( 2P A (0 w k (t) C O S ( 271 f c t)) c( t) M{ ( t) 


k =0 


(4.4) 


where 

P k : Received power from k th channel 
b k (t) : Information bits for k th channel 

w k {t) : Walsh Code for k ,h channel 

In (4.4), k = 1 defines the despread information bits that are meant to be 
recovered from channel. All the other possibilities (k ^ 1) will be considered as despread 
co-channel interference. 

h (0 = v/2if b x (t )w, <( )c( t) cos(27t f c f) c( £ q< 0 
= j2P l b x (1)cos{2Tif c t) 


Y ,(0 


/C-l 


2P k b k (t )w k f )c( 0 cos(2ti ./;./) 


&=0 


v 

AT-l 


c(t)vtf(0 


^(0^,(0 n;(Ocos(27t/ c O 


&=0 

k^\ 


(4.6) 
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Despread Additive White Gaussian Noise, 


n l (t) = n(t)(( t) \\{ t) 


(4.7) 


Even though each VSAT has been separated by different PN sequences, the 
interference between VSATs is unavoidable due to imperfect orthogonality of PN 
functions. Interference increases as the number of active VSATs (denoted as M) is 
increased. Interfering channels from different terminals will possibly contain different 
phase inf ormation (cp ; ) from the channel, which is intended to be recovered in our 
analysis. 


( M K.-l 


Si(0 


EEV24 b lk (t +x,. ) w ik {t +x ( . ) c ; (t +x ; )cos(2ti f c t +cp ( .) 


y /=1 k =0 
M K.-l 




(4.8) 


EE^ b, k (t +x ( ) w ik (t +x,. )w, (t)c,(t+Xf)c(t) cos(2ti f c t + 9.) 


(=1 it=0 


The despread Jammer signal reaches the receiver with no data bit, PN or Walsh 
sequence. The jammer interference is simply a tone located in the middle of the satellite 
receiver bandwidth spectrum. Since the jammer cannot synchronize itself with the 
receiver, a phase difference between the jammer and receiver (<p .) is inevitable. 

cos(2tc f c t+(p.)c(t) w t (t) (4.9) 


2 . The Demodulated Signal, y 2 ( t ) 

A demodulated signal is simply obtained by applying the modulation tone, which 
is coherent detection. The subscript “2” indicates demodulation in our notation. 

y 2 ( 0 = y, O)cos(27t/;0 

= (A (0 + Yi(0 + «i(0 + S i(0 + 7, (t )) cos(27t f c t) (4.10) 

= A(0 +Y2(0 + n 2 (0+G 2 (0+ 7 2 ( t) 

Now, let’s calculate each of the components of the demodulated signal to obtain 
the decision statistic, Y 
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Despread and modulated information bits, I 2 (t ): 

/ 2 = *j2P t k[(t)cos(2jif c t)cos(2nf c t) 

= : ^-^i(0(l+cos(47t/ c O) 

where 

cos 2 (27t/ f t) = ^-(i + cos(4^/ f ,t)). 


Despread and modulated co-channel interference, y 2 (0 : 


Y 2 (0 = £ (t) w k (t ) (0(1 + cos(47t f c t )) 

A=0 


De spread and modulated AWGN 

« 2 (0 = n( t) t{ t) w; (r) cos(27t / c f) 

Despread and demodulated inter-VSAT interference 

M fc 1 I2P 

S 2 (0 = X L,^T~ b ^ t + x i) w ik(t + ^,-)Wi(0 c,(f +T,.)c(0cos(2jt f c t +cp,)cos(27t / c f) 

,=i fc=o 2 

where 

cos(27t f c t +cp ! )cos(27t: //) 

= -^-[cos cp ; -cos(4t / f t)costp ( . +sin(4t / c r)sincp ( .] 

M A',-1 l2p 

S 2 (0 = £ L +T z) +T ,-) (0 <3 (r+x,-M 0 

i=l *=0 ^ 

[coscp,-cos(4t //)coscp, +sin(4t / f t)sincp ; ] 


(4.11) 


(4.12) 


(4.13) 


(4.14) 


44 



Despread and modulated jammer interference 


J 2 (t) = -j2Pj c(t) Mi(^cos(271 / c t + cp^)cos(27t/ c t) 

flP 


c(t ) ^(^[coscfK -cos(47t/ f 7)cos cp ; +sin(47t/ c 7)sincp j .j 


3. The Decision Statistic 

Decision statistic, Y, can be modeled as Gaussian random variable with 
value of Y and variance of c, . The variance of Y will be introduced in next section. 

T 

Y = | y 2 (t ) dt 

0 

Y = J I 2 (t)dt + Jy 2 (0 dt + ^n 2 (t)dt + J<; 2 (t) dt + J J 2 (t)dt 

0 0 _^ 0 _^ 0 _ / 0 _ 

Y n c, J 

where, 

y : Co-channel Interference, 

r K ~ l I2P 

Y = J (Ow k (t)w l (t)(l + cos(4ti f c t )) dt 

o k=0 2 

k *1 


n : Additive White Gauissian Noise 

T 

T] = |/7(t)c( t) H|(T)COS(27t f c t)dt 
0 


(4.15) 


a mean 


(4.16) 


(4.17) 


(4.18) 
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: Inter-VSAT Interference 


i 

q=jq 2 (t)dt 

0 

M K-l j2p 

_ f £ £ VphO (' +T,)»’ 1 (r)c,(r+i,MD 


= 0 


=1 Jk=0 


[coscp, -cos(4t / c ?)costp j + sin(4t f c t)sin(p t ]dt 


M bp f 

= C0S( P, V 0 J b ik (t +X .)w. k (t +T ; .)w i (0 C. (t +X ,M f) (it 


i =1 &=0 


In reference [3] a ik (t +X ( ) and J, (7) are given as 

a ik (t +x,.) = ^ t (f +x,.)w* (t +x i )c i (t +x,.) 
dft) — w ,(t )c ( t) 


(4.19) 


M K.-. 

€=££ 

i =1 fc=0 




coscp 


1 

ja ik (l +x i )d l (/) dl 


i 0 


(4.20) 


/ : Jammer Interference 

r j2P. 

J = J—— -c(t) Wj(f)[coscp y .-cos(47t/ c t)coscp ; . +sin(4t / C t)sincp dt 

Since the multiplication of f c and t was selected as the integer multiplier, 
the integration of cos and sine functions over a period is always zero. 


J = 



T 

coscp 7 |c(t) wft) dt 
0 


(4.21) 


a. Mean Value of Decision Statistic, Y 

t T T T T r ] 

| 1 2 (0 dt + Jy 2 (0 dt + J n 2 (t ) dt + 1<; ft) dt + J J 2 (t) dt \ (4.22) 

0 0 0 0 0 J 
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(4.23) 


Y = ^||/ 2 (r)rfr|+ pA+pA + pA+pA 
= £” |j / 2 (r) ^r| 

_ T I2P 

Y= [ ——- b x (t)(\ + cos(4;t f c t )) dt 
o 2 

JlR 

_ ± V— Lj 
2 

Since the multiplication of f c and t was selected as the integer multiplier, 
the integration of cos and sine functions over a period is always zero. Likewise, 
integration of bft) over a period is ±1. Because the number of ones and minus ones are 
assumed to be equal. 


4. Variance of Decision Statistic, Y 

Hie variance of our decision statistic is the sum of the variances of our co-channel 
interference term and noise plus jammer term, which are assumed to be independent: 

Var {= Var{ y j + Var{n) + Var{q) + Var{j) 

= £{Y ; }+£{« ; }+£{s J } + £p J } (4 ' 24) 

a. Variance of co-channel interference 


Var{y] =£{y 2 ] 


E\ 


IL ( t ) w k ( t ) wft ) (1 + cos(471 f c t )) clt 


IV 


0 *=0 
k*l 


(4.25) 


= 0 


T 

Because J w 1 (t) w k (t) = 0 where k ^ 1. 

0 
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b. Variance ofAWGN 
Var{n ] = £{/z 2 } 

( T T ) 

J/z(0c( t) w\(t)cos(2jif c t)dt^n(k)c(X)w 1 (k)cos(2jif c 'k)d'kl (4.26) 

0 0 J 

( TT t 

JJ 'n(t }i ^i)c(Oc(^)w 1 (Owi(X)cos(27C/ c Ocos(27C f c X)dtdk ? 

We assume that n(t) is a wide-sense stationary white noise process. The 
autocorrelation of the wide-sense stationary white noise process was defined in reference 
[16] as 

E{n(t)n(k)} = ^8(t-X). (4.27) 

so, 


Var{n} = f c 2 (t ) vtf(t)—(l +cos(4 nf c t))dt 

2 o 2 

No T 

4 


(4.28) 


c. Variance of the Inter-VSAT Interference 
Var{q}=E{q 2 } 


M K ,-1 


Var{c ) } = Y J Y, E 


i =1 k—0 


EE 


cos 


1 

( ( ? l )\a ik (t+x i )dft)dt 


Lv 


H ^i-l T T j 

Var{q} = Y,Y, E \ ^f cos2 (9,-) J % +x ; M (0<*f a a (k +x ,)d, (A,) dX [ 


i=l /c=0 


Var{q} 
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M Kr -1 T T 


Var{q}=Ltll fi 2 (t-'k)dtdk 

4 I=1 ^=o o o 

Let’s apply transformation of lim its and variables of the integral 


u = t-X 
v = t + X 

t = ^(u + v) 

X =—(v— n) 
2 V 7 


J* =det 


—> u+v -2t 
—>u —v = —2X 


'at 

ax" 


'1 

r 

3n 

a« 

= det 

2 

2 

3t 

ax 

1 

1 

_ a v 

dv _ 


.2 

2. 


After applying transformation 


M K -—1 p 7] 2T |w| 


'Ms}=££-fJ J pV)4 


cluclv 


i =1 &=0 


m ^-i p ^ 2r :H 

= EE~f j I $ 2 ( u ) dudv 

>= 1 *=o 4 o |„i 


j=i £=o 4 q 


M A/-1 p i 

EEi L JP 2 ( M )( 2T ’- 2w )^ 


P(«) : 


i ^ i ix T 

1-n<7 

rj -r | | C 


0 


otherwise 


M Kj -1 p 


Var{<;} = ££^ | 


r/ . >. 

/y / O A7,. A7"2,.2 \ 


i=l £=0 


, 2 An AV 

1 + —— 


[2T -2u)du 


J 


M Kj -1 p 

=EEy 2 

1=1 *=0 ^ 


„ 2 A 2 n 3 n 2 2An 3 A 2 n 

Tu-Nu~ +-+ 

37 2 


2. .4 


37 47 2 


3o 


EE^ 

1=1 fc=0 ^ 

M Kj -1 p 

= EEy 

i=l fc-0 z ' 


27 2 


-2 y 


3A 6A 


V 3i 7 


for A »1 


(4.29) 
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d. Variance of Jammer Interference 

The derivation of the jammer interference is not much different from the 
inter-VS AT interference, 


Var{j}=E{j 2 } 

[p. 1 (4.30) 

= E j -^-cos 2 ((fv)j j c(f)wft )c Qi)w l (X)dtdX !> 

Recall from (4.19) d l ( t ) = wft)c ( t ). 

p TT 

Var{j} = — 1| E{d l (t)d l (X)}dtdX 

4 0 0 
p T T 

- X) dtd X 

4 0 0 
P T 2r-|u| 

= -[ \ I P (u )J,x dvdu 

-t |«| 
p t 2 r-|«| 

= —J J P (u) dvdu 
4 ° H 

p T 

= — J2(3(w) (T-u) du 

4 o 


P(w) = l-- 


= 1-— where 7 = — 
T N 


P 


Vn 


y ar {y} = _L J 


f 


T —u —Nu + 


Nu 


2 3 


V 


J 


„ M 2 Vw 2 An 3 
7w-h 


iu 

7/ A 
/jv 


37 


_ i 
2 


P, f 7 2 7 2 AT 2 


V 2A- 


- + - 


V7 


J 

'3 A 


2V 2 37V 3 
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_ _J_ 
2 


P . (T 2 


T 2 T 
■ + 


2 3 


N 2N 2N 3 N ■ 


__ j_ 

2 


P ( p 2 


-2 \ 


IN 6N~ 


_ P j T ~ 


AN 


for N »1 


( 4 . 31 ) 


5. Signal-to-Noise Plus Interference Plus Jamming Ratio 


Var\tf\ =Var{n} + Var{c,} + Var{j ) 


N 0 T 

4 



P ik T 2 

6N 


PT 2 

j 

AN 


( 4 . 32 ) 


5/lV/? = 


F 2 

Var{t,} 


p;r 

2 


jv; 

4 


M K-l 

+E£ 


1=1 A:=0 


PJ 2 

6A^ 


+ 


P jP 2 

~Jn 


( 4 . 33 ) 


S/1VR = 


AT;-1 


( 4 . 34 ) 


P,, and fj can be considered equal. The incoming transmissions from the channel 
that is meant to be detected and the interfering channels are assumed to have the same 
power because the size of the footprint on earth is \ery small relative to the geostationary 
satellite altitude. 


SI NR = 


_1_ 

N 0 | M{K -\) | PjT 
2 TP l 3 N 2NPJ 
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PT = E, 


P\T = E„ 


J b _ 


= ,Sy/?(SignaltoJamRatio) 


E 

— = SNR (SignaltoNoiseRatio) 

N A 


SINR 


1 


, M(K- 1) 1 


1 


(4.35) 


2(SiVR) 3iV 


2N(SJR) 


6 . Performance Analysis with Convolutional Encoding 

Proakis states in reference [4] that 


P 2 (d) < Q(yj(d (SINR) R) 


f 

\ 


dR 


1 1 

[i 

2 (SNR) 3 N 2N(SJR) ^ 


f 

\ 


d 


1 1 

ill 

2 R(SNR) 3 NR 2 NR(SJR) 


where R is the code rate. 

The reference [4] defines the Q function as in (4.37). 


(4.36) 


Q(x) = Lerfc(-j=) 


(4.37) 
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And finally, the upper bound for the probability of bit error for the Spread 
Spectrum, DS-CDMA BPSK VSAT system is [4] 

P b < £ M(d) (4-38) 

d =d free 

Where the set of $ d s represents the distance spectmm of the code. 

Therefore, 


n < E 

d~dfr ee 

The probability of bit error equation (4.39) was plotted by using a Matlab® 
simulation. The tone jammer Matlab® simulation calculates the probability of bit error for 
different code rates and for a different number of active channels. The results were 
plotted as probability of bit error, (Pb), versus signal-to-noise ratio ,E b /N 0 . The 
simulation was mn for four different signal power-to-spreading factor times jammer 
power ratio (PJ NP j ). According to the simulation results, as the number of users 

increase, the performance of the VSAT system becomes worse. The performance was 
never acceptable for any traffic intensity and for jammer powers for the coding rates of 
1/8, 1/16 and 1/32 . 




1 


, M(K-l) f 


1 


2R(SNR) 3NR 2 NR (SJR) 


(4.39) 


53 




Figure 19. Performance of VSAT System with the Code Rate of 1/16 
in Low Power Jammer Condition P ] / NP j =15 dB. 

Even in a low power jammer scenario, the codes specified above perform very 
unsatisfactorily. 

Except for those code rates specified cbove, in every other code rate, the system 
can support data communication only for a low number of active users in adverse tone 
jammer conditions ( P l /NP j = 5 dB). N is chosen 256 as previously explained. 

Additionally, lower codes 3/4 and 2/3 did not perform well for severe jamming. Even 
though the best performance was obtained for the code rate of 1/120, there is little 
performance difference between the code rate l/3and 1/120. For a signal-to-noise ratio 
of 15 dB and for a high number of active channels, code 1/120 can perform 0.5 dB better 
than the code rate of 1/3 for an expense of 40 times less bit rate. Therefore code rate 1/3 
is the best solution in a severe tone-jamming situation. 
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Figure 20. Comparison of Code Rates of 1/3 and l/120for P { / A/P = 5 dB 

As the jammer power changes from adverse to upper moderate (P ] / NP i =\0 dB), 

even code rates 3/4 and 2/3 performed well for low loads. However, code 1/3 still had 
the best optimal performance relative to code rate 1/128. Even for the best performance, 

the probability of bit error was insignificant than 7x10 1 . The change in performance for 
different code rates was not too much for heavy loads. This result co nfir ms that as the 
traffic intensity increases, the co-channel interference and inter-VSAT interference are 
the dominant factors on performance. For example, performance differences between the 
code rates of 1/3 and 1/120 diminish as the number of active channels increases in 
Figure 20. 

In low jammer effect conditions ( P l / /VP =15 dB or P l /NP j =20 dB), the 

probability of bit error always remains under 10 7 up to 100 active channels. The best 
performing code rates almost support reliable voice communication in a weak jammer 
scenario. 
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Probability of Bit Error 



Figure 21. Performance of Code Rate 1/3 for P t / NP j = 15 dB . 


Finally, the proposed VS AT system can sustain data communication for all 
jammer-included scenarios via code rate 1/3. Voice communication can only be 
conducted with less active channels in moderate jammer power condition. 
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B. PERFORMANCE ANALYSIS WITH PULSE JAMMER 


Pulse jamming is one of the most effective ECM techniques. In this technique, the 
jammer is activated for a short period of time and paused for the rest of the period. The 
variable duty cycle of the pulse jammer (p ) influences the performance of the targeted 
communication system. The pulse jammer duty cycle (p ) is a function of the ratio of the 
signal power to the total jammer power {PI J) [2]. The lower duty cycle requires much 
bandwidth compared with the higher duty cycle, jamming the wider part of the 
operational bandwidth of the communication system. The jammer can even jam the entire 
bandwidth of the target system by applying the necessary duty cycle. On the other hand, 
since the jammer has l im ited power, the power spectral density ever the unit band would 
be decreased as the jammed band increases. As illustrated in Figure 16 at the beginning 
of Chapter IV, the total area under the power spectral density curve for a jammer will 
always be the same, regardless of the applied duty cycle or the jammed portion of the 
whole band. 


This situation causes the jammer operator to consider a tradeoff between the duty 
cycle and power spectral density pair to force the target communication system to 
perform as poorly as possible. 

As previously mentioned, the jammer is activated for a fraction of the total period. 
The variance of interference that is received by the satellite receiver changes accordingly. 

Jammed bits, r ,, can be modeled as Gaussian random variable with mean -a/eJ volts and 
a variance of, [16] 


(K-l)E c , N 0 , N 


3 V 


+ - 


-+- 


2p 


(4.40) 


where 

Nj = PSD of Jammer 
p = Jammer pulse duty cycle (0 < p < l) 

Unjammed bits, r l2 are also modeled as Gaussian random variable with mean -yj~E~ c as 
well and a variance of, 
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(4.41) 


c2 (K-l)E c ^N { 

3 N e 2 

The combination of two Gauissian variables,! 7 , received by the satellite receiver is 
expected to be Gaussian again. This property was depicted as in (4.42) in reference [4]. 

P 2 (d\i) = plp,>o\ (4.42) 

Because every received bit is not jammed, the equation above can be adapted for 
the noise jammer as, 

P 2 (d\i) = P ||£iu+I>, j2 |>0 (4.43) 

'='+i ) J 

In this expression, we assumed that i of the total d bits a period are jammed and 
(d-i )bits are unjammed. 

1. Mean and Variance of Gaussian Random Variable Y [4] 

i 

M y = Y j m i (4.44) 

/=1 

where m i = 

The jammer will not affect the mean of the Y. 

M, = £ -.Je, + £ -JeT = t,-4 e c 

1 =1 l=i+\ 1 =1 

= ~ds[K 

The variance of r.v. Y is affected by the jammer as explained previously. 
Therefore we must define the variance according to the duty cycle of the jammer. 
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(4.45) 


<q=E°? 

/= 1 

i d 

=E<v ! +E° J 

Z=1 l =i +1 

= /a 2 + (d-i)a 2 


Let’s substitute theo 2 and a 2 found in equations (4.40) and (4.41) respectively. 
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N was described as the total spreading factor. The total spreading factor, N is 
fixed to an integer, 256, which is the multiplication of the spreading factor due to the 
coding and spreading factor due to the PN sequence spreading. In this case, the equation 
(4.46) can be simplified as follows: 

2 iN, d{K-\)E. dN ’ ,. . 

o ;=— L +— - - jL + - 51 (4.47) 

y 2p 31V 2 
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2 . 


Conditional Probability of Error, P 2 (d i ) 


Since we know the exact amount of jammed bits (/) to calculate variance, the 
conditional probability of error P 2 ( d\i ) can be calculated by the following equation [4]. 


P 2 (d\i) = Q 
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(4.48) 


3. First-Event Error Probability, P 2 (d) 

Tlie First-Event Error Probability, P 2 (d) can be found using binomial 
distiibutions. Jamming occurs with the probabihty of p and not being jammed occurs 
with the probabihty of (1-p). In reference [5], PMF of X for the binomial random 
variable X is given as: 
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(4.49) 


(4.50) 


4. Performance Analysis with Convolutional Encoding 

The probability of bit error with convolutional encoding was given in (4.38) in the 
tone jammer performance section. 

Pi,= i p aw 

d~d free 

d ( d \ 

= L PrfE . p i (l-p) d - i P 2 (d\i) (4.51) 
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We can use the same simplification for the term inside the Q function, as we did 
for the tone jamming. 
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The pulse jammer simulation shows similar results to the noise jammer 
simulation. Pulse jammer Matlab® simulation calculates the probability of bit error for 
different code rates, for different numbers of active channels and for different duty cycles 
of the jammer (p). The results were plotted as the probability of bit error (Pb) versus the 

signal-to-noise ratio for each signal-to-jam ratio ( E b /J 0 ) of 5, 10, 15 and 20 dB. 
According to the simulation results, some common patterns were revealed: For example, 
as the code rate increases, system performance becomes more unyielding against the 
variable jammer duty cycle. Deterministically, for the code rate of 1/4, the performance 
of the VSAT system decreases 12 dB as the duty cycle of the jammer decreased to 0.01 
from 0.1 for E h /N 0 = lOdB , E b /J 0 = 15 dB and K = 100 active users (Figure 22). 

However, for the same configuration, the VSAT performance decreases only 0.5 

dB for the code rate of 1/32 (Figure 23). For this reason, choosing a lower code rate 

would be more appropriate to el im inate this capability of the jammer. The expected 
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common pattern for all the results was that the performance of the VS AT system 
becomes worse as the number of users increase due to increasing interference from co¬ 
channels. 

A rate span of 3/4 to rate 1/128 has been studied to observe the effects of a code 
rate over the probability of bit error. Code rates up to 1/4, do not affect the FI severely. 
Figure 24 and Figure 25 illustrate the performance difference between the code rate of 
1/2 and 1/64 for E b /J 0 =15 dB . However, for the rates 1/8, 1/16 and 1/32 the 
performances are not acceptable for all signal-to-jam ratios, similar to the case in tone 
jamming. Then, for the rates 1/64 , 1/128 performance increases again. For instance, for 
the rate of 1/128 the results are slightly better than for the rate 1/2 for p =1 (Figure 26). 
However, rate 1/128 cannot provide data rates as high as rate 1/2 can. The spectrum 
distances algorithm used to obtain the (3 ( , coefficients gives perfect results for the rates of 
1/20, 1/40, 1/60 and 1/120 . Because of these perfect spectral distance coefficients, rate 
1/20 performs better than rate 1/16 . Moreover, rate 1/60 and 1/120 perform better than 
rates 1/64 and 1/128 respectively. 
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Figure 22. Pulse Jammer Effect for Code Rate 1/4 and E h / J 0 = 15 dB 
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Figure 23. Pulse Jammer Effect for Code Rate 1/32 and E b / J 0 = 15 clB 
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Probability of bit error 
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Probability of bit error 



Figure 25. Pulse Jammer Effect for Code Rate 1/64 and E b / J Q = 15 clB 
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Probability of bit error 
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The simulation results also showed that performance is mostly affected by 
intrasystem interference for E b /J 0 values greater than 15 dB. In other words, co-channel 

interference or inter-VSAT interference is more effective than jammer interference for 
this situation. 



Figure 27. Pulse Jammer Effect for Code Rate 1 / 2 and E h / J 0 = 20 clB 


The differences in performance curves between the Figure 24 and Figure 27 
illustrate that changing in performance for the different channel densities decreases as 
the E h / J 0 ratio increases. 

Figure 28 introduces that any signal-to-jam ratio below 5 dB is not enough for a 
reliable communication. In this illustration, one of the best performing code rate, 1/3, 
was used. 

High code rates such as 1/3 or 1/4 guarantees reliable communication for all 
reasonable signal-to-jam conditions (Figure 29). But as the number of active channels 
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Probability of bit error 


increases (over 150), the small duty cycles of the jammer can increase the probability of 
bit error over 10 1 . 
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Figure 29. Pulse Jammer Effect for Code Rate 1 / 3and E h / J 0 =10 dB 


In a jammer free scenario, the rate 1/2 or 1/3 would perform best, enabling all 
users to communicate with high data rates. After the jammer appears in the same 
scenario, the VSAT system would switch to a higher code rate but to a lower data rate 
until communication could be accomplished efficiently. Code 1/20 or 1/40 would be the 
most optimum code rates for moderate jamming. If the jamming is more severe, rate 
1/120 must be used with only data communication. The reason voice communication 
cannot be used is that the proposed VSAT system has a li mited bandwidth. 
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C. SUMMARY 

In this chapter, we modeled the VS AT system by defining the similarities and 
differences between cellular and satellite communication. We have used properties of the 
Walsh function and the autocorrelation function of PN sequences through the whole 
analysis 

Two types of jammer effects were analyzed with Matlab® simulations. We have 
looked for a proposed probability of bit error corresponding to the min im um and 
maximum signal-to-noise ratios, which were introduced in Chapter n. As a result of these 
simulations, the superiority of the pulse jammer over the tone jammer was revealed. In 
other words, the proposed VS AT system is more resistive to the tone jammer than it is to 
the pulse jammer. Considering the hostile side, the worst case for the enemy pulse 
jammer (p = 1) was as effective as the tone jammer. The bwer duty cycles of the pulse 
jammer decreased the performance of our VSAT system dramatically. 

Both jammer capabilities can prevent the platforms from making secure voice 
communication. Except under severe conditions, data communication can be supported 
most of the time. As a countermeasure, decreasing the system bit rate, planning the 
necessary EMCON measures in an operation or decreasing the number of communicating 
channels would even increase the performance more. As previously mentioned in the 
Rain Loss section of Chapter n, increasing the output power of VSAT increases the 
signal-to-noise ratio. Therefore, better performance is obtained. Even though, the overall 
spreading factor was fixed to 256 in this thesis, increasing processing gain by using more 
PN chips per bit increase the durability against enemy jammer effects. 
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V. CONCLUSIONS 


In this thesis, we have presented the tools for a DS-CDMA Spread Spectrum 
communication system using forward error correcting mechanism and analyze the 
performance of proposed system by using those tools. 

Hie main goal in Chapter II was calculating the minimum signal-to-noise ratio 
that our system can support. Initially, we defined the term VSAT and other orbital 
mechanics parameters. Then we decided on the proposed VSAT system requirements. By 
using the physical sizes, distances and power limitations, a link budget analysis was 
made. The worst-case scenario was considered in this analysis assuming a friendly 
platform at the edge of the satellite footprint and the enemy jimmcr at the best place to 
minimize our system performance. The link budget analysis presented the minimum and 
maximum signal-to-noise ratios that can be established at the satellite receiver. The 
difference between the minimum and maximum signal-to-noise ratios has been calculated 
by the rain loss analysis. 

In Chapter IH, we examined the algebraic properties of Walsh functions and PN 
sequences that will help us to figure out how to minimize jammer effects and friendly 
interference on our system. The Walsh functions consisted of orthogonal sequences, 
which support full coverage between sequences. We also presented the extended 
orthogonality that ensures the orthogonality between every channel in a VSAT. The PN 
sequences were used to spread out our frequency spectrum, reducing the power spectral 
density, and minimizing the jammer effects. Afterwards, we introduced the convolutional 
forward error correction coding, which lets the system perform better with a low signal- 
to-noise ratio. 

The performance analyses were introduced in Chapter IV. The performance 
analyses were considered both for tone jammer and pulse jammer. First of all, we 
developed a signal-to noise plus interference ratio and the probability of bit error for our 
DS-CDMA satellite system by using Gaussian approximation. We used Gaussian 
approximation to define median and variances of random variable Y . For both jammer 
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cases, co-channel interference, inter-VSAT interference, Additive White Gaussian Noise 
and Jammer Interference were considered as variances of Y . We introduced Forward 
Error Correction in the form of convolutional encoding with soft-decision decoding into 
our uplink model and developed an upper bound on the probability of bit error [3]. 

We used Matlab® to mn the simulations. The simulation results for both jammer 
types revealed that the proposed VS AT uplink system could accommodate the data or 
voice communication within reasonable lim its. The performance plots were given in 
probability of bit error versus signal-to-noise ratio for the different jammer powers and 
the number of active channels. We used the obtained minimum signal-to-noise ratio in 
Chapter II to decide whether system performance was sufficient or insufficient. The 
following results are valid for both pulse and tone jammer. 

• As the number of active channels increases, the probability of bit error 
increases. 

• As the jammer power increases, the probability of bit error increases. 

• There was no direct correlation between the code rate and the system 
performance. 

• The performances of code rates of 1/8, 1/16 and 1/32 were unacceptable. 

• The code rates 1/3,1/120 and 1/128 were best performing code rates. 

• The code rate 1/3 introduces most efficient performance and bit rate pair. 

Additionally, some results were subject to specific jammer types. 

• The worst case for enemy pulse jammer (p =1) affected VS AT uplink as 
badly as the tone jammer. 

• As the duty cycle of pulse jammer decreases, the probability of bit error 
increases. 

• As the pulse jammer power decreases, the effect of changing duty cycle 
decreases. In other words, for the high E b / J 0 values effects of tone 
jammer and pulse jammer are almost the same. 
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Finally, we must recall that we never gave up considering the worst-case scenario 
in our analysis. In the real word we would expect better performance from our VS AT 
system. Also, changing antenna size, increasing the VSAT power or decreasing bit rate 
can always empower the VSAT system against unpredictable and variable jammer 
capabilities. 
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